Novel Fermi Liquid of 2D Polar Molecules 
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We study Fermi liquid properties of a weakly interacting 2D gas of single-component fermionic 
polar molecules with dipole moments d oriented perpendicularly to the plane of their translational 
motion. This geometry allows the minimization of inelastic losses due to chemical reactions for 
reactive molecules and, at the same time, provides a possibility of a clear description of many-body 
(beyond mean field) effects. The long-range character of the dipole-dipole repulsive interaction 
between the molecules, which scales as 1/r^ at large distances r, makes the problem drastically dif- 
ferent from the well-known problem of the two-species Fermi gas with repulsive contact interspecies 
interaction. We solve the low-energy scattering problem and develop a many-body perturbation 
theory beyond the mean field. The theory relies on the presence of a small parameter kpr^, where 
kp is the Fermi momentum, and r« = md^ /h^ is the dipole-dipole length, with m being the molecule 
mass. We obtain thermodynamic quantities as a series of expansion up to the second order in kpr-t 
and argue that many-body corrections to the ground-state energy can be identified in experiments 
with ultracold molecules, like it has been recently done for ultracold fermionic atoms. Moreover, 
we show that only many-body effects provide the existence of zero sound and calculate the sound 
velocity. 



I. INTRODUCTION 

The recent breakthrough in creating ultracold diatomic 
polar molecules in the ground ro- vibrational state [H-Q 
and cooling them towards quantum degeneracy [l| has 
opened fascinating pro spects for the observation of novel 
quantum phases [5l422l|. A serious problem in this direc- 
tion is related to ultracold chemical reactions, such as 
KRb+KRb=>K2-l-Rb2 observed in the JILA experiments 
with KRb molecules [l^, [l^l, which places severe lim- 
itations on the achievable density in three-dimensional 
samples. In order to suppress chemical reactions and per- 
form evaporative cooling, it has been proposed to induce 
a strong dipole-dipole repulsion between the molecules 
by confining them to a (quasi)two-dimensional (2D) ge- 
ometry and orienting their dipole moments (by a strong 
electric field) perpendicularly to the plane of the 2D 
translational motion [2^, [2^ . The suppression of chem- 
ical reactions by nearly two orders of magnitude in the 
quasi2D geometry has been demonstrated in the recent 
JILA experiment [l^]. At the same time, not all polar 
molecules of alkali atoms, on which experimental efforts 
are presently focused, may undergo these chemical reac- 
tions [1^. In particular, they are energetically unfavor- 
able for RbCs bosonic molecules obtained in Innsbruck 
or for NaK and KCs molecules which are now being 
actively studied by several experimental groups (see, e.g. 
[29}). It is thus expected that future experimental studies 
of many-body physics will deal with non-reactive polar 
molecules or with molecules strongly confined to the 2D 
regime. 

Therefore, the 2D system of fermionic polar molecules 
attracts a great deal of interest, in particular when they 
are in the same internal state. Various aspects have been 



discussed regarding this system in literature, in partic- 
ular the emergence and beyond mean field description 
of the topological Px + ipy phase for microwave-dressed 
polar molecules ll|, Tsj. in.terla.yer sup erfluids in bilayer 
and multilayer systems [l^, llJ, [l^ ISOf. the emergence of 
density- wave phases for tilted dipoles [l6l[T7l . [20l - |22j . The 
case of superfluid pairing for tilted dipoles in the quasi2D 
geometry beyond the simple BCS approach has been dis- 
cussed in Ref . [12] ■ The Fermi liquid behavior of this sys- 
tem has been addressed by using the Fourier transform 
of the dipole-dipole interaction potential [13, [H, IsTI - 
Is^ l and then employing various types of mean field ap- 
proaches, such as the Hartree-Fock approximation [l3] or 
variational approaches [l^, [11] . It should be noted, how- 
ever, that the short-range physics can become important 
for the interaction between such polar molecules, since in 
combination with the long-range behavior it introduces 
a peculiar momentum dependence of the scattering am- 
plitude [18]. 

On the other hand, there is a subtle question of many- 
body (beyond mean field) effects in the Fermi liquid be- 
havior of 2D polar molecules, and it can be examined in 
ultracold molecule experiments. For the case of atomic 
fermions, a milestone in this direction is the recent result 
at ENS, where the experiment demonstrated the many- 
body correction to the ground state energy of a short- 
range interacting two-species fermionic dilute gas (35l . [36| . 
This correction was originally calculated by Huang, Lee, 
and Yang [st', 'ss'] by using a rather tedious procedure. 
Later, it was found by Abrikosov and Khalatnikov [39] in 
an elegant way based on the Landau Fermi liquid theory 

In this paper, we study a weakly interacting 2D gas 
of fermionic polar molecules which are all in the same 
internal state. It is assumed that each molecule has an 
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average dipole moment d which is perpendicular to the 
plane of the translational motion, so that the molecule- 
molecule interaction at large separations r is 

U{r) ^ \ ^ (1) 

where r* = jtkP' jfi' is the characteristic dipole-dipole 
distance, and m is the molecule mass. The value of d 
depends on the external electric field. At ultralow tem- 
peratures that are much smaller than the Fermi energy, 
characteristic momenta of particles are of the order of the 
Fermi momentum kp^ and the criterion of the weakly in- 
teracting regime is: 

kpr^ < 1. (2) 

As a consequence, the Fermi liquid properties of this 
system, such as the ground state energy, compressibility, 
effective mass, can be written as a series of expansion in 
the small parameter kpr^,. We obtain explicit expressions 
of these quantities up to the second order in kpr^,, which 
requires us to reveal the role of the short-range physics 
in the scattering properties and develop a theory beyond 
the mean field. Our analysis shows that only many-body 
(beyond mean field) effects provide the existence of un- 
damped zero sound in the coUisionless regime. 

The paper is organized as follows. In Section II we an- 
alyze the low-energy 2D scattering of the polar molecules 
due to the dipole-dipole interaction. We obtain the scat- 
tering amplitude for all scattering channels with odd or- 
bital angular momenta. The leading part of the ampli- 
tude comes from the so-called anomalous scattering, that 
is the scattering related to the interaction between par- 
ticles at distances of the order of their de Broglie wave- 
length. This part of the amplitude corresponds to the 
first Born approximation and, due to the long-range 1 /r^ 
character of the dipole-dipole interaction, it is propor- 
tional to the relative momentum k of colliding particles 
for any orbital angular momentum /. We then take into 
account the second Born correction, which gives a con- 
tribution proportional to k'^. For the p-wave scattering 
channel it is necessary to include the short-range contri- 
bution, which together with the second Born correction 
leads to the term behaving as k'^ In k. In Section III, af- 
ter reviewing the Landau Fermi liquid theory for 2D sys- 
tems, we specify two-body (mean field) and many-body 
(beyond mean field) contributions to the ground state 
energy for 2D fermionic polar molecules in the weakly 
interacting regime. We then calculate the interaction 
function of quasiparticles on the Fermi surface and, fol- 
lowing the idea of Abrikosov-Khalatnikov [39] , obtain the 
compressibility, ground state energy, and effective mass 
of quasiparticles. In Section IV we calculate the zero 
sound velocity and stress that the many-body contribu- 
tion to the interaction function of quasiparticles is neces- 
sary for finding the undamped zero sound. We conclude 
in Section V, emphasizing that the 2D gas of fermionic 
polar molecules represents a novel Fermi liquid, which 



is promising for revealing many-body effects. Moreover, 
we show that with present facilities it is feasible to ob- 
tain this system in both coUisionless and hydrodynamic 
regimes. 

II. LOW-ENERGY SCATTERING OF 
FERMIONIC POLAR MOLECULES IN 2D 

A. General relations 

We first discuss low-energy two-body scattering of 
identical fermionic polar molecules undergoing the 2D 
translational motion and interacting with each other at 
large separations via the potential U{r) ([T]). The term 
low-energy means that their momenta satisfy the inequal- 
ity kr^ <C 1. In order to develop many-body theory for 
a weakly interacting gas of such molecules, we need to 
know the off-shell scattering amplitude defined as 

/(k',k)= y exp(-ik'r)C/(r)4(r)d2r, (3) 

where ipklj) is the true wavefunction of the relative mo- 
tion with momentum k. It is governed by the Schrodinger 
equation 

A + C/(r)Uk(r) = ^k(r). (4) 

\ m J m 

For |k'| — |k| we have the on-shell amplitude which enters 
an asymptotic expression for V'k(r) at r — oo l4l|: 



V'k(r) = exp(ikr) - -^^f^^^ exp(ifcr), (5) 

with <p being the scattering angle, i.e. the angle between 
the vectors k' and k. 

The wavefunction '0k (r) can be represented as a sum of 
partial waves ipi{k,r) corresponding to the motion with 
a given value of the orbital angular momentum /: 

oo 

■!/'k(r) = ^ ■0;(fc, r)i'exp(i;(/3). (6) 
Using the relation 

oo 

exp(jkr) = ^ i'^Ji{kr)exp[il{ifk-'fr)], (7) 

^— — oo 

where J; is the Bessel function, and (fk and (fr are the 
angles of the vectors k and r with respect to the quan- 
tization axis. Eqs. ^ and ([7]) allow one to express the 
scattering amplitude as a sum of partial-wave contribu- 
tions: 

oo 

/(k',k) = eM^l^)fiik\k), (8) 

/ — — oo 
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with the off-shell /-wave amplitude given by 



B. Low-energy p-wave scattering 



fi{k\k)^ Ji{k'r)U{r)iji{k,r)2T:rdr. (9) 
Jo 

Similar relations can be written for the on-shell scattering 
amplitude: 



(10) 



fi{k)= Ji{k'r)U{r)i:i{k,r)2TTrdr. (11) 

The asymptotic form of the wavefunction of the /-wave 
relative motion at r — )■ oo may be represented as 



cos(fcr — 7r/4 + 5i{k)) 



kr 



(12) 



where 5i{k) is the scattering phase shift. This is obvious 
because in the absence of scattering the /-wave part of the 
plane wave exp(ikr) at r — )■ oo is {kr)~^/'^ cos(fcr — 7r/4). 
Comparing Eq. (fT2|) with the /-wave part of Eq. ([5]) we 
obtain a relation between the partial on-shell amplitude 
and the phase shift: 



4/1^ iAn5i{k) 
ml — ita,Ti5i{k) 



(13) 



Note that away from resonances the scattering phase shift 
is small in the low- momentum limit fcr* ^ 1. 

For the solution of the scattering problem it is more 
convenient to normalize the wavefunction of the radial 
relative motion with orbital angular momentum / in such 
a way that it is real and for r — > cxd one has: 

Mk,r)^ [Ji{kr) - i&n6i{k)Ni{kr)] 

cxcos(fcr-/7r/2-7r/4 + (5i(fc)), (14) 

where Ni is the Neumann function. One checks straight- 
forwardly that 



^i{k,r) 



ipi{k,r) 



I — i taiiSi{k) 



Using this relation the off-shell scattering amplitude 
can be represented as 



(15) 



where fi{k',k) is real and follows from Eq. (|9]) with 
ipi{k,r) replaced by il)i{k^r). Setting k' — k we then 
obtain the related on-shell scattering amplitude: 



4r 

fi{k,k) = fi{k)^ tan5,(fc). (16) 

m 



As we will see, the slow decay of the potential 
U(r) at sufficiently large distances makes the scattering 
drastically different from that of short-range interacting 
atoms. For identical fermionic polar molecules, only the 
scattering with odd orbital angular momenta / is possi- 
ble. For finding the amplitude of the p-wave scattering 
in the ultracold limit, kr^, <C 1, we employ the method 
developed in Ref. [18] and used there for the scattering 
potential containing an attractive dipole-dipole tail. 
We divide the range of distances into two parts: r < 7'o 
and r > ro, where ro is in the interval r* <C ro <C fc~^. In 
region I where r < tq , the p-wave relative motion of two 
particles is governed by the Schrodinger equation with 
zero kinetic energy: 



1 dil)i 
r dr 



0. 



(17) 



At distances where the potential U{r) already acquires 
the form ([T]), the solution of Eq. ([T7)) can be expressed in 
terms of growing and decaying Bessel functions: 



V'/(r) oc 




(18) 



The constant A is determined by the behavior of U{r) at 
short distances where Eq. ([1} is no longer valid. If the 
interaction potential U{r) has the form ([1]) up to very 
short distances, then A = 0, so that for r -> equation 
()18|) gives an exponentially decaying wavefunction. 

It should be noted here that for the quasi2D regime ob- 
tained by a tight confinement of the translational motion 
in one direction, we can encounter the situation where 

^ 'oj with /g being the confinement length. However, 
we may always select rg ^ lo if the condition fc/g ^ 1 is 
satisfied. Therefore, our results for the 2D p-wave scat- 
tering obtained below in this section remain applicable 
for the quasi2D regime. The character of the relative mo- 
tion of particles at distances r ^ /g is only contained in 
the value of the coefficient A, and the extra requirement 
is the inequality /c/g ^ 1. 

At large distances, r > rg, the relative motion is prac- 
tically free and the potential U{r) can be considered as 
perturbation. To zero order, the relative wavefunction is 
given by 



Ji(fcr) - ia.n5i{k)Ni{kr), 



(19) 



where the phase shift 5i{k) is due to the interaction be- 
tween particles in region I. Equalizing the logarithmic 
derivatives of ipi (r) and ilif) at r = rg we obtain: 



tan 5 1 — ^ 



Tvk'^ror^: 
8 



l + H ( 2C---2A+ln- 
ro \ 2 To 



(20) 



with C — 0.5772 being the Euler constant. 
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We now include perturbatively the contribution to the 
p-wave scattering phase shift from distance r > tq. In 
this region, to first order in U{r), the relative wavefunc- 
tion is given by 

= rpfhr)- / G{r,r')U{r')rpf}{r')2nr'dr', (21) 



where the Green function for the free p-wave motion 
obeys the radial equation: 

m \ dr^ r dr J ' 27rr 

For the normalization of the relative wavefunction chosen 
in Eq. (IT4l) . we have: 



G(r,r') 



4h' 



^l)f}{r')Ni{kr), r>r' 
Jjfj'{r)Ni{kr'). r<r' 



(22) 



Substituting this Green function into Eq. 1^1^ and taking 
the limit r oo, for the first order contribution to the 
phase shift we have: 



[^fi\r)fU{r)2TTrdr. 



Using Eqs. (fT9|) and ([20l) we then obtain 



(23) 



tan5(^'(fc) = -^-^^ (^2A+2C+ln^-^-^ . (24) 

To second order in U{r), we have the relative wave- 
function: 

/>oo 

^f^{r)=%l;f}{r)+ I G{r,r')U{r')2Trr' dr' 



X / G{r',r")U{r")^pfj>{r")2Trr"dr". (25) 



Taking the limit r — ;> oo in this equation we see that 
including the second order contribution, the scattering 
phase shift becomes: 

tan5i(fc)=tan5(i)(fe)-^ / {rfUir)2nrdr 

8h Jro 

POO 

X / Ni{kr')U{r')'iljfj> {r')2Trr'dr' . (26) 

J r 

As we are not interested in terms that are proportional 
to or higher powers of fc, we may omit the term 
tani5/(fc)A^i(fcr) in the expression for tpfj'ir). Then the 
integration over dr' leads to: 

tan<5i(fc) = Uii S['\k) - ^^^^ r ^^dx 

2 Jkro X 



-x{No{x)J2{x) - Ni[x)Ji{x)) 



1 



']:Nq{x)Ji{x) + \Ni{x)J2{x) 

2 D TTX 



(27) 



For the first four terms in the square brackets, we may 
put the lower limit of integration equal to zero and use 
the following relations: 



dr 1 

/ Jl[x)N,{x)'^ = -^, 
Jo X An 



J!{x)J2{x)No{x)- 



dx 1 



J!{x)Nq{x) 



dx 



167r' 



r°° dr 1 

Jl{x).h{x)N,{x)-^ = - — . 

For the last term in the square brackets we have: 
r ,2, ,dx 1 C ln2 1 , , 

We then obtain: 

tanSiik) ^t£,nS['\k)-^^^^ 
2kr^ nk'^r'^ 



(28) 
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■ In ^fcr. 



C - In 2 + In fcro 
(29) 



where: 



f = exp ( 3C-ln2 - - 2A 



(30) 



Using Eqs. and ([^ we represent the on-shell p- 
wave scattering amplitude fi{k) in the form: 



/i(fc) = /«(fc) + /r(fc). 



(2), 



with 



and 



/«(fc) = |^fc.. 
3 m 



fi'\k)^^{kr,)^\nikr,. 
2m 



(31) 



(32) 



(33) 



The leading term is f[^\k) oc k. It appears to first order 
in U{r) and comes from the scattering at distances r ~ 
1/k. This term can be called "anomalous scattering" 
term (see [4lj). The term /p-* (fc) cx fc^ In^fcr, comes from 
both large distances ~ 1/k and short distances. Note 
that the behavior of the wavefunction at short distances 
where U{r) is no longer given by Eq. ([1]), is contained in 
Eq. ([2^ only through the coefficient <^ under logarithm. 



C. Scattering with |i| > 1 

The presence of strong anomalous p-wave scattering, 
i.e. the scattering from interparticle distances ~ 1/A:, 



5 



originates from the slow \/r^ decay of the potential U{r) 
at large r. The strong anomalous scattering is also 
present for partial waves with higher /. In this section 
we follow the same method as in the case of the p-wave 
scattering and calculate the amplitude of the Z-wave scat- 
tering with \l\ > 1. For simplicity we consider positive I, 
having in mind that the scattering amplitude and phase 
shift depend only on \l\. 

To zero order in U(r), the wavefunction of the /-wave 
relative motion at large distances r > ro is written as: 

= [Mkr) - ti,ii6i^j)ik)Ni{kr)] , (34) 

where (5;(/)(fc) is the /-wave scattering phase shift coming 
from the interaction at distances r < vq. We then match 
'ip^jj-^{k^r) at r = ro with the short-distance wavefunc- 
tion ipi(i)(j') which follows from the Schrodinger equation 
for the /-wave relative motion in the potential U{r) at 
k = 0. This immediately gives a relation: 



tan(5j(7)(fc) 



kJ'iikro) - wiJijkro) 
kN[(kro) ~ wiNiikroY 



(35) 



where the momentum-independent quantity wi is the 
logarithmic derivative of il^i{i){r) aX r — vq. Since we 
have the inequality kr^ <C 1, the arguments of the 
Bessel functions in Eq. ([55)) are small and they reduce 
to Ji{x) ~ cc' , Ni{x) ^ This leads to tanSnr){k) ~ 
(fcro)^'. Thus, the phase shift coming from the inter- 
action at short distances is of the order of (/crp)^'. As 
we confine ourselves to second order in fc, we may put 
tan(5;(/)(fc) = for the scattering with |/| > 1. 

Then, like for the p-wave scattering, we calculate the 
contribution to the phase shift from distances r > ro 
by considering the potential U{r) as perturbation. To 
first and second order in J7(r), at r > ro we have similar 
expressions as Eq. (|23| , ([25)1 for the relative wavefunction 
of the /-wave motion. Following the same method as in 
the case of the p-wave scattering and retaining only the 
terms up to /c^, for the first order phase shift we have: 



tan 5'i'^\k) 



4^2 / [^l"h)ir)?U{r)2nrdr 



nkr^ 



/ Jl{x) — dx 

J kro 



4/2 



(36) 



The second order phase shift is: 



tan5p^(fc) 



™' /(o) ^^^2 



%lu){rYU{r)2-Krdr 



X / iVi(fcr')[/(r')V'i°A)(r')27rr'dr' 



{^kr^Y p Jf{x) 



dx 



kro 



-2 ^^2/' (37) 



and we may put the lower limit of integration equal to 



zero. For the integral over dy, we obtain : 
Ni{y)My) 



yZ 



-dy 



^ :Ji{x)Ni^i{x) + ^j^-^Ji+i{x)Ni{x) 



2/(2/- 1) 

+ ^^[Ni-i{x)Ji+i{x)-Mx)Ni{x)] - ^. (38) 
Then, using the relations: 



X3 4/(/2-l)' 



oo t2 



Ni^i{x)Ji+i{x)dx = 



4/(/-fl)7r' 



:lli^Ni{x)dx = - ^ 



4/27r' 



-A^Ji{x)Ni-i{x)dx = 



8/2(/-M)7r' 



:^Ji+iix)Ni{x)dx^ 

X 



8/2(/ + l)7r' 



we find the following result for the second order phase 
shift: 



tan5p)(fc) = M^ 



1 



' " 8 /(/2- l)(4/2- 1)' 

So, the total phase shift is given by 

tan Si{k) = tan Sl^\k) + tan Sl^\k) 



(39) 



37r(fcr, 



4/2-1 8/(/2-l)(4/2-l)' 



(40) 



Then, according to Eq. ([TB)) the on-shell scattering am- 
plitude fi{k) is 



Mk) = ff\k) + ff\k), (41) 



where 



-(1) , 8?i fcr* 1 



frik) 



fT\k) = -^{kr^Y 



4/2-1' 



1 



2m ' |/|(/2- l)(4/2- 1)' 



(42) 



(43) 



Note that Eqs. (HU and do not contain short-range 
contributions as those are proportional to A:2|'l and can 
be omitted for 1/1 > 1. 
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D. First order Born approximation and the leading 
part of the scattering amplitude 



As we already said above, in the low-momentum limit 
for both 1^1 = 1 and |^| > 1 the leading part of the on-shell 

scattering amplitude fi{k) is fi^\k) and it is contained 
in the first order contribution from distances r > ro- For 
|Z| > 1 it is given by Eq. dHJ and follows from Eq. 

with V';^//) — Ji{kr). In the case of |Z| = 1 this leading 
part is given by Eq. (j32p and follows from the integral 
term of Eq. (I23p in which one keeps only Ji(fcr) in the 
expression for il)fj{r). This means that f^^^k) actually 
follows from the first order Born approximation. 

The off-shell scattering amplitude can also be repre- 
sented as fiik', k) = fl^\k', k) + fl'^\k\ k), and the lead- 
ing contribution J^'^^ {k' , k) follows from the first Born ap- 
proximation. It is given by Eq. Q in which one should 
replace ipi{k,r) by Jiikr): 

poo 

fl^\k\k)= i Ji{kr)Ji{kr')U{r)2Tirdr. (44) 

"'0 

Note that it is not important that we put zero for the 
lower limit of the integration, since this can only give a 
correction which behaves as fc^ or a higher power of k. 
Then, putting U{r) — h^r^jmr^ in Eq. (HU, we obtain: 



(45) 



where F is the hypergeometric function. The result of 
Eq. (1451) corresponds to k < k' , and for k > k' one should 
interchange k and k' . 

For identical fermions the full scattering amplitude 
contains only partial amplitudes with odd I. Since the 
scattered waves with relative momenta k' and — k' corre- 
spond to interchanging the identical fermions, the scat- 
tering amplitude can be written as (see, e.g. |4]|): 

/(k',k) -/(k',k)-/(-k',k). (46) 

Then, according to equation (|10p one can write: 

(47) 



/(k',k) = 2^/,(fc',fc)exp(zM- 



I odd 

In the first Born approximation there is no difference 
between fi{k', k) and Ji{k' , k) because tan5i(fc) in the de- 
nominator of Eq. is proportional to k and can be dis- 
regarded. Therefore, one may use fl^\k' ,k) of Eq. ps)) 
for fi{k',k) in Eq. (|T7|) . One can represent /(k',k) in a 
different form recalling that in the first Born approxima- 
tion we have: 



/(k',k) = y l/(r)exp[i(k-k')r]d^ 



(48) 



Performing the integration in this equation, with U{r) 
given by Eq. ([!]), and using Eq. PSI) wc obtain: 



/(k',k) = 



-{|k + k'|-|k-k'|}. (49) 



Equation (j49p is also obtained by a direct summation over 
odd / in Eq. (gT]), with fi{k', k) following from Eq. (gS]). 



III. THERMODYNAMICS OF A WEAKLY 
INTERACTING 2D GAS OF FERMIONIC 
POLAR MOLECULES AT T = 

A. General relations of Fermi liquid theory 

Identical fermionic polar molecules undergoing a two- 
dimensional translational motion and repulsively inter- 
acting with each other via the potential ([T]) represent a 
2D Fermi liquid. General relations of the Landau Fermi 
hquid theory remain similar to those in 3D (see, e.g. 
|40|). The number of "dressed" particles, or quasipar- 
ticles, is the same as the total number of particles iV, 
and the (quasi)particle Fermi momentum is 



47riV 



(50) 



where S is the surface area. At T = the momentum 
distribution of free quasiparticles is the step function 



n(k) = e{kF - fc). 



(51) 



i.e. ?T.(k) = 1 for k < kp and zero otherwise. The chemical 
potential is equal to the boundary energy at the Fermi 
circle, fjL = ep = e{kp). 

The quasiparticle energy e(k) is a variational deriva- 
tive of the total energy with respect to the distribution 
function n(k). Due to the interaction between quasipar- 
ticles, the deviation 6n of this distribution from the step 
function (|5ip results in the change of the quasiparticle 
energy: 



Mk)= I i^(k,k')<5n(k')^ 



(52) 



The interaction function of quasiparticles i^(k, k') is thus 
the second variational derivative of the total energy with 
regard to n(k). The quantity 6n(k) is significantly dif- 
ferent from zero only near the Fermi surface, so that one 
may put \i = kpn and k' = kpii' in the arguments of 
F in Eq. ([5^ . where n and n' are unit vectors in the 
directions of k and k' . The quasiparticle energy near the 
Fermi surface can be written as: 



eik) = ep + nvFik-kp)+ I F(k,k')Mk')7^' 



2- (53) 



The quantity vp = de{\f.) / Tidk\}.^]^p is the Fermi velocity, 
and the effective mass of a quasiparticle is defined as 
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m* = hkp/vF- It can be obtained from the relation (see 



B. Two-body and many-body contributions to the 
ground state energy 



1 

m 



1 

TO* 



1 



(27r?i)2 



2tt 



F{d) COS 9d9, (54) 



where is the angle between the vectors n and n', and 
F{0) = F{kFn,kFn'). _ 
The compressibihty k at T = is given by [40] : 



(55) 



The chemical potential is /x = ep, and the variation of 
II due to a change in the number of particles can be ex- 
pressed as 



F{kFn,k')Sn{k') 



d^k' 



dk, 



■Ski 



(56) 



The quantity (5n(k') is appreciably different from zero 
only when k' is near the Fermi surface, so that we can 
replace the interaction function F by its value on the 
Fermi surface. Then the first term of Eq. ([55]) becomes 



m 



de 



,5n(k'~ 



d'^k' 



SN 



■27r J "■^'"'(27r)2 2TrS 
The second term of Eq. ([55)) reduces to 



F{e)de. 



dkp m 
We thus have (see [40|): 



-okp = —OKp — 



TO 



S 



(57) 



dN 



1 f^"" , , 2TTh^ 
2ttS Jo "J 



2TTh^ 



1 



iS 2tiS 



{1 -cos 9) F{0)d0. (58) 



Equation ([58]) shows that the knowledge of the interac- 
tion function of quasiparticles on the Fermi surface, F(9), 
allows one to calculate dfi/dN and, hence, the chemical 
potential = dE/dN and the ground state energy E. 
This elegant way of finding the ground state energy has 
been proposed by Abrikosov and Khalatnikov [331 ■ 1^ 
was implemented in Ref. [39j for a two-component 3D 
Fermi gas with a weak repulsive contact (short-range) 
interspecies interaction. 

We develop a theory beyond the mean field for cal- 
culating the interaction function of quasiparticles for a 
single-component 2D gas of fermionic polar molecules in 
the weakly interacting regime. We obtain the ground 
state energy as a series of expansion in the small param- 
eter kpr^t: and confine ourselves to the second order. In 
this sense our work represents a sort of Lee-Huang- Yang 
[37I Issj and Abrikosov-Khalatnikov [sl] calculation for 
this dipolar system. As we will see, the long-range char- 
acter of the dipole-dipole interaction makes the result 
quite different from that in the case of short-range inter- 
actions. 



We first write down the expression for the kinetic en- 
ergy and specify two-body (mean field) and many-body 
(beyond mean field) contributions to the interaction en- 
ergy. The Hamiltonian of the system reads: 



k 



k°k+2^ J2 t^(q)aki+q«L-q«k2aki,(59) 



ki,k2,q 



where aj^. and Ok are creation and annihilation opera- 
tors of fermionic polar molecules, and U{q) is the Fourier 
transform of the interaction potential U{r): 



U{q) = J d^vU{r)e-''' 



(60) 



The first term of Eq. ([55)1 represents the kinetic energy 
and it gives the main contribution to the total energy 
E of the system. This term has only diagonal matrix 
elements, and using the momentum distribution (|5ip at 
T = we have: 



E, 



kin 



s 



n^k'^ 2Trkdk 
2m (27r)2 



h^k% 

16to 



(61) 



The interaction between the fermionic molecules is de- 
scribed by the second term in Eq. ()59|) and compared 
to the kinetic energy it provides a correction to the to- 
tal energy E. The first order correction is given by the 
diagonal matrix element of the interaction term of the 
Hamiltonian: 

^^'^^^ E t^(q)(«l+qaL-q«k.ak,) 



ki,k2,q 



^ J2 [C/(0)-C/(k2-ki)]nk,nk2. (62) 



25* 



ki,k2 



The second order correction to the energy of the state | j) 
of a non-interacting system can be expressed as: 



_g(2) _ VjmVmj 



m^3 



(63) 



where the summation is over eigenstates |to) of the non- 
interacting system, and Vjm is the non-diagonal matrix 
element. In our case the symbol j corresponds to the 
ground state and the symbol to to excited states. The 
non-diagonal matrix element is 



ki,k2,q 



c/(q)«L- 



^k2-q"k2aki 



j). (64) 



This matrix element corresponds to the scattering of two 
particles from the initial state ki , k2 to an intermediate 
state k'j^, k2, and the matrix element Vmj describes the 
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reversed process in which the two particles return from 
the intermediate to initial state. Taking into account 
the momentum conservation law ki + k2 = k'^ + kj the 
quantity Vj„iV„ij = |V,mp is given by 



1 



(2^)^ 



X |C/(k'i-ki)-C/(k^-ki)|', 



(65) 



and the second order correction to the ground state en- 
ergy takes the form: 



£;(2) 



1 



(25)2 ^ 



|C/(k'i-ki)-[/(k'2-ki) 



"-kl?^k2(l - '^k'J(l - «ki) 



(66) 



From Eq. (pS)) we see that the second order correc- 
tion diverges because of the term proportional to riki 'T-k2 j 
which is divergent at large k'l . This artificial divergence is 
eliminated by expressing the energy correction in terms of 
a real physical quantity, the scattering amplitude. The 
relation between the Fourier component of the interac- 
tion potential and the off-shell scattering amplitude is 
given by [4lj : 



U{k' -k")fik",k) 



(67) 



where i?k = h?k^ /ra and i?k" — h^k"^/m are relative 
collision energies. Obviously, we have: iJk — ^/k" = 
n'^iki^ + ka^ - kf - kf )/2m, with ki, ka {k'l, k^') 
being the momenta of colliding particles in the initial 
(intermediate) state, as the relative momenta are given 
by k = (ki - k2)/2, k" = (k'l' - k^')/2. We thus can 
write: 



;7(k'-k)-/(k',k)- 



2m ^ J7(k'-k")/(k",k) 
Wsj-f\^i + kl-kf-k'^-tO 



(68) 



Then, putting k' = k we have 

2to ^ C/(k- k")/(k",k) 



C/(0) = /(k,k) 



h^s j-fkf + kl- k'l 



L-//2 
^2 



and setting k' = — k we obtain 



U(k k)-f(kk) 2m ^ [/(-k-k")/(k",k) 

^1 

Using these relations the first order correction (p^ takes 
the form: 



= ^ E [/(k,k)-/(-k,k)]nk,nk. 

ki,k2 

[;7(k-k')-[/(-k-k')]/(k',k) 

'^Slj::^^,f{kl + ^-kf-k',^-^0)/2m 



?iki»T-k2- (69) 



The quantity [U{k - k') - U{-k - k')] in the second 
term of Eq. ()69p . being expanded in circular harmonics 
eicp{ilip) contains terms with odd I. Therefore, partial 
amplitudes with even / in the expansion of the multiple 
/(k',k) vanish after the integration over cPk'. Hence, 
this amplitude can be replaced by [/(k', k) — /(k', — k)]/2. 
As we are interested only in the terms that behave them- 
selves as ~ fc or ~ fc2^ the amplitudes in the second term 
of Eq. (|69l) are the ones that follow from the first Born 
approximation and are proportional to k. Therefore, we 
may put [C/(k- k') - C/(-k- k')] = [/(k, k') - /(-k, k')] 
and /(k',k) = /*(k, k'). Then the first order correction 
takes the form: 

^^'^ ^ ^ E [/(k, k) - /(-k, k)] nk,nk2 - 



E 



ki,k2,ki 



ki.k2 

|/(k',k)~/(k',~k)|2 

, h'{kl + kl-kf-k'i-i0)/2j 



-nki^k^. (70) 



Using the expansion of the full scattering amplitude 
in terms of partial amplitudes as given by Eq. (|T7|) we 
represent the first order correction as 

i?W = 4 W/.(fc)nk,r.k.-^EE 



5 ^ ^ 

ki ,k2 / odd 



ki ,k2 / odd 



<fk' 



ffik) 



-rikinkj. (71) 



(27r)2 n'{kl+kl-k'}-k'^~i{))/27 

The contribution of the pole in the integration over (Pk' 
in the second term of Eq. ([TT]) gives imff{k)/4,h? for each 
term in the sum over ki, k2, and /, and we may use here 
the amplitude fl^\k). In the first term of Eq. ((7T|) we 
should use fi{k) = fi^\k) + fi^\k). However, we may 

replace //^'' by //^'' because the account of tand{k) in 
the denominator of Eq. ([T5|) leads to k^ terms and terms 
containing higher powers of k. For the amplitude fj^^k), 
we use the expression: 



which assumes a small scattering phase shift. The second 
term of this expression, being substituted into the first 
line of Eq. ([7T]) . exactly cancels the contribution of the 
pole in the second term of (|7T|) . Thus, we may use the 
amplitude // in the first term of equation (I71|) and take 
the principal value of the integral in the second term. The 
resulting expression for the first order correction reads: 

^'^^4E/>k."«^^-(^ 

X V 2r^|/(k,k') /( k,k')|2 

k.iki n\kl+kl-kf-k'i) 

where /(k) = oddU^)- ^ 

The second order correction (j66p can also be expressed 

in terms of the scattering amplitude by using Eg. (1671) . 



f(i) 



(1) 



F(i)i2, 
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Replacing U(ki - k[) = C/(k - k') and - ki) = 
?7(— k — k') by /(k',k) and /(— k, k'), respectively, we 
have: 



i5(2) 



1 



s [ 



l/(k^k)-/(k^-k)|^ 



ki,k2,ki 

X f^kl?^k2(l - '^k'j)(l - ^k^) 



(73) 



Note that the divergent term proportional to nkink2 
in Eq. ([73|) and the (divergent) second term of Eq. ((72|) 
exactly cancel each other, and the sum of the first and 
second order corrections can be represented as E^^^ + 
^(2) ^^(1)+ ^(2) where 



?^kl"k2^ 



(74) 



ki,k2 



and 



l/(k',k)-/(k',-k) 



(2^)' k.t^^ki ^ ^'(k? + ki - kf - kf )/2m 

X nkinkj(l - 7iki)(l - "-ki) - !]}■ (75) 

The term E^^^ originates from the two-body contribu- 
tions to the interaction energy and can be quoted as the 
mean field term. The term E^'^^ is the many-body con- 
tribution, which is beyond mean field. 

It is worth noting that the term proportional to the 
product of four occupation numbers vanishes because its 
numerator is symmetrical and the denominator is anti- 
symmetrical with respect to an interchange of ki , k2 and 
k']^,k2. The terms containing a product of three occu- 
pation numbers, ?T.k2 '^k'j and rtkifT-k2"ki are equal to 
each other because the denominator is symmetrical with 
respect to an interchange of k[ and k2 . We thus reduce 
Eq. dZSl) to 



^(2)=_ 



_j_ ^ 2m|/(k^k)-/(k^~k)| 

2^1 tTk' h\kl+kl-kf-k'i) 



-nki"-k2"-ki-(76) 



Equations ([7^ and ([75]) allow a direct calculation of the 
ground state energy. With respect to the mean field 
term E^^^ this is done in Appendix [Xj However, a di- 
rect calculation of the many-body correction E^^^^ is even 
a more tedious task than in the case of two-component 
fermions with a contact interaction. We therefore turn 
to the Abrikosov-Khalatnikov idea of calculating the 
ground state energy (and other thermodynamic quan- 
tities) through the interaction function of quasiparticles 
on the Fermi surface. 



C. Interaction function of quasiparticles 

The interaction function of quasiparticles F{k, k') is 
the second variational derivative of the total energy with 



respect to the distribution n^. The kinetic energy of 
our system is linear in rik (see Eq. ([5^ ). and the second 
variational derivative is related to the variation of the 
interaction energy E. We have [40| : 



(77) 



k,k' 



where E = E^^^ + E^'^\ and the quantities and M"^^ 
are given by equations (|74p and (|76p . On the Fermi sur- 
face we should put |k| = |k'| = kp, so that the interaction 
function will depend only on the angle 9 between k and 
k'. Hereinafter it will be denoted as F{9). 

The contribution F'-^\9) = 2SdE'^^^/6ni,Snk' is given 

by 



F'-^\9) = 2/ 



|k-k' 



/ odd 



= 2j2filkF\sin-\), (78) 



where /; — //^^ -t- and the amplitudes //^^ and 

//^^ follow from Eqs. ([321) and ([Ml) at |/| = 1, and from 
Eqs. (gH), (021) at |/| > 1. We thus may write equation 



frik) = 

for any odd /, and 



/ln(?fcn); 



\ |i|(i2-l)(4/2-l) ' 1^1 > 1 



with ^ from Eq. ([3D]). Making a summation over all odd 
I we obtain: 



I odd 



/ odd 



2^ 

ln(efcr,)-Y^ + 31n2 



(79) 



,(80) 



Putting k = kp \ sin(6'/2)| and substituting the results of 
equations and dBU]) into Eq. ^ we find: 



F'- >(9) = |sm-| H (/cFf*) 

m 2 m 

9 25 
ln|^r,fcFsin-| - — + 31n2 



(81) 



The many-body correction (l76l) we represent as E^^^^ 
|ki-kip 



e[^^ where 



p(2)^ 8(7r?tr^)^ ^ 

1 mS'2 ^ k^ + ki-kf-k, 

ki.k2,k; 12 12 

^(2)^ SjTrhr^r J2 |ki-k'iHk2-k'i| 



Ta^kinka^k;, (82) 



mS^ ^ kf+ki-kf-kf"'^^"'^^"''^' 

ki,k2,ki 12 12 
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and we used Eqs ((46l) and (|49| for the scattering ampli- 
tudes. The contribution to the interaction function from 
i?^^'' is calculated in Appendix 151 and it reads: 



The contribution from E^^^ is calculated in Appendix [Cl 
It is given by 



9 /4 



37r+27rsin Inl tan - 

2 l3 ' 2 



(84) 



71" , , ^'i , , ^'i , , , ^'ix ^/.s 4arcsin|sin|| - 27r 
^ln2 + --7rln|sin-|+41n|cos-|-41n(l + |sin-|)+g(0) + ^—^ 



TT — 2 arcsin I sin — I 
' 2' 



COS -In- 



|sin(0/2)| 



|sin(0/2)| 



2|sin-| 



(85) 



where 



Q{9) = j 2 sin^ iy9 In f sin + ^in^ ^ + cos^ ^ sin^ (/3 ) dc/j, (86) 



so that 



dg{e) 9 |sin(0/2)| 1 arcsin|cos(0/2)| / 9 „, 

~ir " 2 2 " sin(0/2) cos(^/2) + | cos(^/2) | y'''' 2 " 2 J ' ^^^^ 



We thus have F{9) = F<-^^9)+F^^\9)+F^^\9), where 
\ ^ follow from Eqs. §1^,, and (P)). 

This allows us to proceed with the calculation of ther- 
modynamic quantities. 



D. Compressibility, ground state energy, and 
effective mass 

We first calculate the compressibility at T = 0. On the 
basis of Eq. ([58]) we obtain: 



^^l _ 27rr 1 
'on ~ mS ^ 2^ 



iO) F^^^9) + F^^\9) + F^^\9) 



d9 



mS 3mS 



mS 



irmS 



(fcFn)^(30-8G-h21C(3)), (88) 



where G — 0.915966 is the Catalan constant, and C(3) = 
1.20206 is the Riemann zeta functi on. Cal culating coef- 
ficients and recalling that kp = y^AnN/S we represent 
the inverse compressibility following from Eq. (j55p in a 
compact form: 

2m 5 \ 37r 4 / 



where we obtain the coefficient Ci — 2.16exp(— 2A) by 
using Eq. (150]) for the coefficient ^ which depends on 
the short-range behavior through the constant A (see 
Eq. (IT8l) V For the chemical potential and ground state 
energy we obtain: 
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^ " + ^^^""^ (^ln[4efcf - - j + — -^(fc^r,)2 - __(fc^^,)2(30 _ 8G + 21C(3)) 

= ^ (l + ^fci^r. + ^{kFr. f ln(C2fcFr0) ■ (90) 

= ^ + ii'^^* + J(fc.^*)^ln(C3fc..,)) , (91) 



with numerical coefficients ^2 = 1.68exp(— 2A) and C3 — searching for the variations Sn in the form: 
1.43exp(— 2^). Note that the first term in the second 
hne of Eq. ([55]) and the first terms in the first fines of 
Eqs. (|5D|) and (PT|) represent the contributions of the ki- 
netic energy, the second and third terms correspond to 
the contributions of the mean field part of the interaction 
energy, and the last two terms are the contributions of 
the many-body effects. 

The effective mass is calculated in a similar way by 
using Eq. (|54)) : 



Sn{q, r, t) — 5{eq — eF)v{n) exp j(kr — ujt) 

and using Eq. ([5^ . from the kinetic equation in the col- 
lisionless regime: 



d5n 
~dt 



dSn driq dSeq 
dr dq fidr 



= 0, 



1 



1 



1 



m* ni (27r?i)2 
1 



(f(i) (0) ^ {o)+f!^^'^ (9)) cos ede 



m 
1 

m 



1 -I- ^kpr^ + ^ikpr*)'^ hi(C4fcFf*) 
Stt 4 



(92) 



where the numerical coefficient ^4 = 0.65 exp(— 2^1). 
Note that if the potential U{r) has the dipole-dipole form 
(H)) up to very short distances, we have to put A = in 
the expressions for the coefficients (i, (2, Ca? C4- Con- 
sidering the quasi2D regime, this will be the case for r, 
greatly exceeding the length of the sample in the tightly 
confined direction, Iq. Then, as one can see from equa- 
tions dHU), dll]), and dll]), the terms proportional 
to (kpr^,)'^ are always negative in the considered limit 
fci?r, <^ 1. These terms may become significant for 
fcpr* > 0.3. 



IV. ZERO SOUND 

In the collisionless regime of the Fermi liquid at very 
low temperatures, where the frequency of variations of 
the momentum distribution function greatly exceeds the 
relaxation rate of quasiparticles, one has zero sound 
waves. For these waves, variations Sn{q, r, t) of the mo- 
mentum distribution are related to deformations of the 
Fermi surface, which remains a sharp boundary between 
filled and empty quasiparticle states. At T — >• the equi- 
librium distribution riq is the step function (jSip . so that 
driq/dq — —nSlq—kp) = —hw5{tq — eF), where v = vpn, 
with n being a unit vector in the direction of q. Then, 



one obtains an integral equation for the function v{n) 
representing displacements of the Fermi surface in the 
direction of n 

k f 

(uj — vpn ■ k)u{n) = - — • k / F(kpn, kpn')i'(n')dn' . 
{2nYTi J 

Introducing the velocity of zero sound Uq — ui/k and 
dividing both sides of this equation by vpk we have: 

777* cosfl /■^'^ - 
[s ~ cos e)v{B) - ^^-^ / F{d ~ B')v{B')de', (93) 

where s = u^/vp^ and 9, 9' are the angles between k and 
n, n', so that — 6*' is the angle between n and n'. The 
dependence of the interaction function of quasiparticles 
F = +F^^^+F^'^'> on {9-9') follows from Eqs. ([8T|) . 
([Ml) , and dSS]) in which one has to replace 6 by {6 - 9'). 

The solution of equation ([55)) gives the function i'{9) 
and the velocity of zero sound uq, and in principle one 
may obtain several types of solutions. It is important to 
emphasize that undamped zero sound requires the con- 
dition s > 1, i.e. the sound velocity should exceed the 
Fermi velocity [40l|. We will discuss this issue below. 

For solving Eq. (j93p we represent the interaction func- 
tion _F as a sum of the part proportional to kpr^. and the 
part proportional to {kpr^)'^. As follows from Eqs. (|8ip . 
I, and (|55|). we have: 



F[9-9')^^kpr, 



sm ■ 



-J^(kpr:)^<^{9-9')m 
m 



where the function ^{9 — 9') is given by the sum of three 
terms. The first one is the term in the second line of 
Eq. (j8ip . the second term is the expression in the square 
brackets in Eq. (j84p . the third term is the one in curly 
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brackets in Eq. (|85l) . and we should replace 9hy {9—0') in 
all these terms. It is important that the function ^{6—6') 
does not have singularities and $(0) — <i>(±27r) = 2-k. 
Using Eq. the integral equation (^51) is reduced to 
the form: 



cose)v{e) 



(3 cos ( 
2m* 



9') 



d0' 



■ COSf 



e / i'ie')^{e-e')d9', 



where /? = (m* /Trm)kpr^ <C 1. 

We now represent the function 1/(9) as 



HO) 



Cp cos J 

p=0 



(95) 



(96) 



Then, integrating over d9' in Eq. (|95|) . multiplying both 
sides of this equation by cos j9 and integrating over d9, 
we obtain a system of linear equations for the coefficients 
Cj. We write this system for the coefhcients r/j — Cj{l 

- 1/4)), so that Cj = r;j(l + 
{p — 1/4 — /3)~"^. The system reads: 

[s - 1)(1 + /3f/o)77o + [r;o - \rii] + PU^m = y^o; (97) 



PUj), where Uj — 



4)(l + /3C/i)?7i + [77i 



1 



- 

?7o- 2'?2]+/3C/im=Y$i;(98) 



(s-l)(l+/3C/,)?7,+[r;,-i(ry,_i+7?^ij+/3C/,77,=^$,; j>2, (99) 
where 



■ / COS 6* cos jOdO I Cp 
Jo ./o p^o 



cosp6i'$(0 - 9')de(lQQ) 



with Cj = 1 for i > 1 and Cq = 1/2, and we put m* — m 
in the terms proportional to 

In the weakly interacting regime the velocity of zero 
sound is close to the Fermi velocity and, hence, we have 
(s - 1) < 1 (see, e.g. ^). Since /3 < 1, we first 
find coefficients r]j omitting the terms proportional to /3 
and in Eqs. For j > 1 equation ^ then 

becomes: 



(» 



1),, - 



2 dp 

and searching for s > 1 we may write 



0, 



exp{-v/2(s-l)j}; j » 1. 



(101) 



If j <C l/y/s — T, then we may also omit the terms pro- 
portional to (s — 1) in the system of linear equations for 
Vj ([97)1 - (|99|) . The system then takes the form: 



0; 



0; 



Without loss of generality we may put ryo = 1/2. This 
immediately gives rjj — 1 for j > 1, which is consistent 
with Eq. (jlOip at j <^ l/^/s — 1. We thus have the zero 
order solution: 



% = 1/2; 

^j = l; 1<j<1/n/F^. 



(102) 



In order to find the coefficients rjj taking into account 
the terms linear in /3, we consider j such that j3Uj ^ 
P/P ^ (s— 1), i.e. j <^ [3/{s — 1). Then we may omit 
the terms proportional to (s — 1) in equations (P7|) -(|M 1) . 
Omitting also the terms proportional to 0^ this system 
of equations becomes: 



m - + /3t^o»7o = 0; 
Vi-ila- i^m + f^UiTji = 0; 



(103) 
(104) 



V, - 2 Nj-i + ^J+i] + PUjTjj = 0; J > 2. (105) 

Putting again 779 — 1/2 the solution of these equations 
reads: 



?7i = 1 + /3C/o 
'7j 



l + /3jC/o + 2/?^(j-p);7p; j>2. 



Confining ourselves to terms linear in /3 we put Up = 
^/{p^ — 1/4) and, hence, Uq = —4. Then, using the 
relation 



1 



_ 4(j - 1) 
^ p2 „ 1/4 2j _ 1 ' 

which is valid for j > 2, we obtain: 



771 = 1 - 4/3; 



(106) 



For j > •\//3/(s — 1) we should include the terms pro- 
portional to (s — 1) in Eq. (|105p . This leads to the 
solution in the form of the decaying Bessel function: 
V] - v/2(s - l)/7rifyjy^(Vs - Ij), which for small /3 

is practically equivalent to Eq. (jlOip . 

Wc now make a summation of equations (in71) - (|99p from 
J = to J = ^ l/ys— 1. The summation of the sec- 
ond terms of these equations gives •y/(s — l)/2, whereas 
the contribution of the terms proportional to (s — 1) is 
much smaller and will be omitted. The sums X^iLo ^i'Hj 
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and X]j=o converge at j ^ l/Vs — 1, and the up- and equation (I107P becomes: 
per limit of summation in these terms can be formahy 
replaced by infinity. We thus obtain a relation: 



'-^+Y^I3^^U,-P^J2^,=0. (107) 

3=0 3=0 



Confining ourselves to contributions up to /3^, in the 
second term on the left hand side of Eq. (jl07l) we use 
coefficients rjj given by Eqs. (|106p . and write Uj = 
l/(j^ - 1/4) + - 1/4)^. In the expressions for $j 

we use Co — 1/2 and Cp = 1 for p > 1. We then have: 

oo oo ^ 

J2 PVjUj = -2/3+J2 {8+51-2^2- 2^3} -(108) 

j=o 1=1^ ' 



The contribution linear in /3 vanishes because 
E^iVO'^ - 1/4) = 2. The quantities 81,82, 



and 5*3 are given by 

00 -j^ 



TT^ -8; 



1 

1 NT^ P 



^2 Ej-2_i/4Ep2_i/4 E(j-2_i/4)(^-+i/2); 



^(j-l/2)(j2-l/4)^ 



so that 



^2 + ^3 = E 



2j = 



^ (j2 - 1/4)2 2 ■ 

We thus see that the contribution quadratic in /3 also van- 
ishes because the term in the curly brackets in Eq. (jlOSp 
is exactly equal to zero. Hence, we have J^JLo I^VjUj — 
up to terms proportional to . 

The sum in the third term on the left hand side of 
Eq. (fTOT)) . after putting Cq = 1/2 and Cp = 1 for p > 1 
in the relations for $ j , reduces to 

^ l>j= — ^ / cos cos jOdO 

0=0 



47r — .,0 



j = -oo 

00 p2-IT 



X 



/ cosp9'^{e-e')de'. (io9) 

p=-oo'^0 

For 9 in the interval < 6* < 27r we have a relation: 

00 

cosje = n[S{9) + S{9-2n)], 

j=-oc 

which transforms Eq. (|109p to 

00 

Y = ^[2$(0) + $(27r) + 4'(-27r)] = 27r2, (110) 

3=0 



2„2 



0. 



This gives s = 1 -|-2(/37r)^, and recalling that = kpr^,/Tr 
(we put m* = to) we obtain for the velocity of zero sound: 



ua — vf[1 + 2{kpr^,Y 



(111) 



Note that in contrast to the 3D two-species Fermi 
gas with a weak repulsive contact interaction (scattering 
length a), where the correction {uq — vp) exponentially 
depends on kpa, for our 2D dipolar gas we obtained a 
power law dependence. This is a consequence of dimen- 
sionality of the system. 

It is important that confining ourselves to only the 
leading part of the interaction function F, which is 
proportional to kpr,, and is given by the first term of 
Eq. (|8ip . we do not obtain undamped zero sound (s > 1) 
[42| . This corresponds to omitting the terms fP^j/2 in 
equations (l97l) - (|99|) and is consistent with numerical cal- 
culations [22l. Only the many-body corrections to the 
interaction function of quasiparticles, given by equations 
and ([55)) . provide non-zero positive values of $(0) 
and $(±27r), thus leading to a positive value of {uq — vp). 
One then sees that many-body effects are crucial for the 
propagation of zero sound. 

In principle, we could obtain the result of Eq. (|llip 
in a simpler way, similar to that used for the two- 
species Fermi gas with a weak repulsive interaction (see, 
e.g. HO])- Representing the function v{9) as v{9) = 
cos 9v{9) /(s — cos 9) we transform Eq. (|93p to the form: 



(27r?i)2 







'^nO-0>iO')cosO'^,_ (112) 
s — cos 9' 



Since s is close to unity, it looks reasonable to assume 
that the main contribution to the integral in Eq. (|112p 
comes from 9' close to zero and to 27r. Using the fact 
that F{9) = F{2tt - 9) we then obtain: 



m*F{9)i>{0) 
4^ 



s- 1 



(113) 



We now take the limit 9^0 and substitute F{0) = 
{^TrTi^ /m){kpr^Y as follows froms Eqs. dH]), and 
([55)1 . Putting TO* = m we then obtain s = 1 + 2(fci?r*)'* 
and arrive at Eq. (|llip . 

Note, however, that for very small 9 or 9 very close 
to 27r the dependence F{9) is very steep. For — > the 
leading part of the interaction function, which is linear in 
kpr^, vanishes, and only the quadratic part contributes 
to F{Q). Therefore, strictly speaking the employed pro- 
cedure of calculating the integral in Eq. (|112l) is ques- 
tionable for very small 9. This prompted us to make the 
analysis based on representing v{9) in the form (196]) and 
on solving the system of linear equations (|n7| - (|99p . 

Equation (jll2p is useful for understanding why un- 
damped zero sound requires the condition s > 1 so that 
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uq > vp- For s < 1 there is a pole in the integrand of 
Eq. (|112p , which introduces an imaginary part of the in- 
tegral. As a result, the zero sound frequency uj will also 
have an imaginary part at real momenta k, which means 
the presence of damping (see, e.g. [iojV 

We could also consider an odd function i^(0), namely 
such that iy{2TT - 9) ^ ~v{e) and vifi) = v{2i:) = 0. In 
this case, however, we do not obtain an undamped zero 
sound. 



V. CONCLUDING REMARKS 

We have shown that (single-component) fermionic po- 
lar molecules in two dimensions constitute a novel Fermi 
liquid, where many-body effects play an important role. 
For dipoles oriented perpendicularly to the plane of trans- 
lational motion, the many-body effects provide signifi- 
cant corrections to thermodynamic functions. Revealing 
these effects is one of the interesting goals of up-coming 
experimental studies. The investigation of the full ther- 
modynamics of 2D polar molecules, including many-body 
effects, can rely on the in-situ imaging technique as it 
has been done for two-component atomic Fermi gases 
[H,!!!]. This method can also be extended to 2D systems 
for studying thermodynamic quantities [i^ |3l ■ Direct 
imaging of a 3D pancake-shaped dipolar molecular sys- 
tem has been recently demonstrated at JILA [45|. For 
2D polar molecules discussed in our paper, according to 
equations (IM1) - (PT|) . the contribution of many-body cor- 
rections proportional to (fcj^r,)^ can be on the level of 
10% or 20% for fc^f* close to 0.5. Thus, finding many- 
body effects in their thermodynamic properties looks fea- 
sible. 

It is even more important that the many-body effects 
are responsible for the propagation of zero sound waves 
in the coUisionless regime of the 2D Fermi liquid of po- 
lar molecules with dipoles perpendicular to the plane of 
translational motion. This is shown in Section IV of our 
paper, whereas mean-field calculations do not find un- 
damped zero sound [l^l- Both coUisionless and hydro- 
dynamic regimes are achievable in on-going experiments. 
This is seen from the dimensional estimate of the relax- 
ation rate of quasiparticles. At temperatures T <C 
the relaxation of a non-equilibrium distribution of quasi- 
particles occurs due to binary collisions of quasiparticles 
with energies in a narrow interval near the Fermi surface. 
The width of this interval is ^ T and, hence, the relax- 
ation rate contains a small factor {T/epY (see, e.g. (40|'). 
Then, using the Fermi Golden rule we may write the in- 
verse relaxation time as r^^ ~ {g1^^/h){m/TT?)n{T/eY, 

where n is the 2D particle density, the quantity ^ m/fi^ 
represents the density of states on the Fermi surface, and 
the quantity ffe// is the effective interaction strength. 
Confining ourselves to the leading part of this quantity, 
from Eqs. (17^ and (17^ we have geff ^ h^kFr^/m. We 



thus obtain: 

{kpr,f[ — \ . 114 

Note that as ep « h'^kp/2m « 27r?i^n/m, for considered 
temperatures T <^ ep the relaxation time r is density 
independent. Excitations with frequencies a; <C 1/t are 
in the hydrodynamic regime, where on the length scale 
smaller than the excitation wavelength and on the time 
scale smaller than l/w the system reaches a local equi- 
librium. On the other hand, excitations with frequen- 
cies oj ^ l/r are in the coUisionless regime. Assuming 
T ^ lOnK, for KRb molecules characterized by the dipole 
moment d ~ 0.25 D in the electric field of 5kV/cm as ob- 
taind in the JILA experiments, we find r on the level of 
tens of milliseconds. The required condition T <C e_F is 
satisfied for ^ 70 nK, which corresponds to n > 2 • 10^ 
cm~^. In such conditions excitations with frequencies of 
the order of a few Hertz or lower will be in the hydrody- 
namic regime, and excitations with larger frequencies in 
the coUisionless regime. 

The velocity of zero sound is practically equal to the 
Fermi velocity vp = fikp/m*. This is clearly seen from 
Eq. (jllip omitting a small correction proportional to 
(fc^r*)'*. Then, using Eq. (1^^ for the effective mass and 
retaining only corrections up to the first order in kpr^,, 
we have: 

hkp f A , \ , , 

Uo^ — {l + —kpn . (115) 

In the hydrodynamic regime the sound velocity is: 

In ^kp ( % \ 

V m oN m \ 6n J 

where we used Eq. (|88p for dfi/dN and retained correc- 
tions up to the first order in kpr,,. The hydrodynamic 
velocity u is slightly larger than the velocity of zero sound 
uq, and the difference is proportional to the interaction 
strength. This is in sharp contrast with the 3D two- 
component Fermi gas, where uq ~ vp > u !v vp/y/3. 

We thus see that it is not easy to distinguish be- 
tween the hydrodynamic and coUisionless regimes from 
the measurement of the sound velocity. A promising way 
to do so can be the observation of damping of driven ex- 
citations, which in the hydrodynamic regime is expected 
to be slower. Another way is to achieve the values of kpr^, 
approaching unity and still discriminate between uo and 
u in the measurement of the sound velocity. For example, 
in the case of dipoles perpendicular to the plane of their 
translational motion the two velocities are different from 
each other by about 20% at kpr^, ~ 0.5. These values 
of kpr^ are possible if the 2D gas of dipoles still satis- 
fies the Pomeranchuk criteria of stability. These criteria 
require that the energy of the ground state correspond- 
ing to the occupation of all quasiparticle states inside 
the Fermi sphere, remains the minimum energy under an 
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arbitrarily small deformation of the Fermi sphere. The 
generalization of the Pomeranchuk stability criteria to 
the case of the 2D single-component Fermi liquid with 
dipoles perpendicular to the plane of their translational 
motion reads: 



2ir 



F{9) COS j9 d9 > 0, (117) 



and this inequality should be satisfied for any integer j. 
As has been found in Ref . 22^ , the Pomeranchuk stabil- 
ity criteria (|117p are satisfied for kpr^ approaching unity 
from below if the interaction function of quasiparticles 
contains only the first term of Eq. ([5T|) . which is the 
leading mean field term. We have checked that the sit- 
uation with the Pomeranchuk stability does not change 
when we include the full expression for the interaction 
function, F{9) = F^^\9) + f[^\9) + F^^\9), following 
from Eqs. (gT]), dUl), and (|HS]). Thus, achieving kpr* 
approaching unity looks feasible. For KRb molecules 
with the (oriented) dipole moment of 0.25 D the value 
kpr^ K, 0.5 requires densities n~ 2 ■ 10^ cm~^. 

Finally, we would like to emphasize once more that our 
results are applicable equally well for the quasi2D regime, 
where the dipole-dipole length r* is of the order of or 
smaller than the confinement length /q = (?i/™'^o)^^^j 
with being the frequency of the tight confinement. 
The behavior at distances r ^ Iq \a contained in the co- 
efficient A defined in Eq. p^ . Therefore, the results for 
the velocity of zero sound which is independent of are 
universal in the sense that they remain unchanged when 
going from ^ Zq to < Iq. The only requirement is 
the inequality kpl^ <^ 1. It is, however, instructive to 
examine the ratio r*/Zo that can be obtained in exper- 
iments with ultracold polar molecules. Already in the 
JILA experiments using the tight confinement of KRb 
molecules with frequency wo ~ 30 kHz and achieving the 
average dipole moment d ~ 0.25 D in electric fields of 5 
kV/cm, we have ~ 100 nm and Iq ~ 50 nm so that 
r*/Zo — 2. A decrease of the confinement frequency to 5 
kHz and a simultaneous decrease of the dipole moment 
by a factor of 2 leads to r*/Zo ~ 0.2. On the other hand, 
for d close to 0.5 (which is feasible to obtain for other 
molecules) one can make the ratio r^t/Zo close to 10 at 
the same confinement length. 
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Appendix A: Direct calculation of the first order 
contribution to the interaction energy 



For directly calculating the first order (mean field) con- 
tribution to the interaction energy i?^^) ([7i|) . we represent 
it as = ^ 



E''2'' where 



(1) 



E, 



(1) 



/(I) 

/(2) 



;(1) 

|ki 



2 

Iki-kal 



d^kid^k2 

l2^ 



(Al) 
(A2) 



and the amplitudes /^^^ and /^^^ are given by Eqs. ([7^ 
and (|80|) . respectively. In the calculation of the integrals 



for E- 



(1) 



and ijj^-* we turn to the variables x = (ki 



k2)/2fcF and y = (ki +k2)/2kF, so that d'^kid^k^ = 
Snkpd'^xd'^ydif, where (p is the angle between the vectors 
X and y, and the integration over dip should be performed 
from to 2tt. The distribution functions and n^^ ^-^e 
the step functions (j5ip . The integration over dki and dk2 
from to kp corresponds to the integration over dy from 
to yoix, (p) = —x\ cosiy9| -|- ■y/r^^"x^"siii2^ and over dx 
from to 1. Using Eq. ([79)) we reduce Eq. (IA1[) to 



(1) 



where 



/i = 



El 



dp I x^dx 



sn^k% 

— ^ kpr^Ii, 



ydy = ^ 



2tt 



dip 



(A3) 



^dx 



X [1 - 2| cos ip\\Jl — x'^ sin^ ip + x'^{cos^ ip — sin^ ip)]. 

The last term of the second line vanishes, and the inte- 
gration of the first two terms over dip and dx gives: 



h = 



TT — 2x\/ 1 — 



2 arcsin x]=—. 

45 



Then Eq. ^M^ yields 



E 



(1) 



8S h^k% 



457r-^ 



128 n^k'^ 



S 45 



(A4) 



which exactly coincides with the second term of the first 
line of Eq. (PT|) . 

Using Eq. (|80p the contribution E^'' takes the form: 



25 

ln(efcj^r*)- — + 31n2 



2iT'^m 

where the integrals I2 and are given by 



/2+/3kA5) 



27r pi 

dip 

Jo 



x^dx 



yo{x.,v) 



ydy 



1 /■^'^ /-^ 

— / dp / X dx 

"2 Jo Jo 



X [1 — 2| cosip\x\J I — x^ sin^ ip + x^{cos^ ip — sin^ ip) 

^ pi ^ 

— — x^\2Tr — 4a;vl — a;^ — 4arcsinxl(ia; = — , 
2 Jo J 32' 



16 



and 

I3. — I dip I X In xdx 



ydy 



X a; 



lna::[l — 2| costplxy 1 — sin^ + x^(cos^ — sin^ ip 

— j x^ lnx[27r — 4a;-\/l — — 4 arcsinxjda; 
2 Jo 




Mi--'' 



Substituting the calculated I2 and I3 into Eq. (jAsp we 
obtain: 



^ 647rm 



ln(4^fcFf*) 



23 
12 



23 

ln(4efcFr*)- — 



(A6) 



This exactly reproduces the third term of the first line of 
Eq. dH]). 



FIG. 1. (color online). Left: The integration area for Ji (in 
blue). The distance between the points Oi and P is RoiP = 
k/2fcF, Rpo2 = k'/2kF, and RoiN = ki/2A:F. Right: The 
integration area for I2 and I3 (in red). The distance between 
the points O2 and A'' is Ro2iv = ki/2fcF, RoiP = k/2fcF, 
and R02P ~ k'/2fcF- 



The integration region is shown in Fig|TJ where the dis- 
tance between the points Oi and O2 is II01O2 = 
The distance between the points Oi and is RqiAt = 
ki/2A:F, and RAro2 = k2/2A:F, so that Ron — x/2. The 
quantity |x| changes from to li{p) where 

6 6 
liif) + cos^ - — 2li{ip) cos - cos ip = 1, 



Appendix B: Calculation of the interaction function 



(2) 



~ (2") 

The interaction function F-^ ' is the second variational 
derivative of the many-body contribution to the interac- 
tion energy, E\ ' ([5^ . with respect to the momentum 
distribution function. It can be expressed as 



Ff)(k,k') 



2h' 



{kFr,f[h+i2+h), (Bl) 



where 



|ki|<fcj!. '^F 

d^ki 



|k-k' 



|ki|<fej=. >^F 



d^ki Iki - k 



/|2 



^K.^^kf + k^-k'^-ki'^'^^*-'''-''- 



and the presence of the Kronecker symbols dq reflects 
the momentum conservation law. On the Fermi surface 
we put |k| = |k'| = kp and denote the angle between k 
and k' as 9. Due to the symmetry property: J^(k, k') — 
i^(k',k) we have F{6) — F(27r — 9) and may consider 9 
in the interval from to vr. 

In order to calculate the integral Ii , we use the quan- 
tities s = (k -f k')/2/eF and m = (k — k')/2fcF and turn 
to the variable x = (ki — k2)/2fcF = (2ki — s)/2kp. For 
given vectors k and k', the vectors s and m are fixed and 
|s| = cos(6'/2), |m| = sin(0/2). The integral can then be 
rewritten as: 



.d^x. 



and li{(p) ■ li{ip + it) = sin^(0/2), with ip being an angle 
between m and x. In the polar coordinates the integral 
Ii takes the form: 



/l = / dp^ 

Jo Jo 



-1 + 2sin^ 



1 



2sin2(0/2)-a;2 

and after a straightforward integration we obtain 

- f 9 9 
Ii ~ n \ 2 sin^ - In I tan - 1 — 1 



xdx. 



(B5) 



In the integral I2, using the variable y = (ki -f k2)/2fcF 
we observe that it changes from to hi'p) where 

9 9 
l^ip) + sin^ - — 2l2{p>) sin - cos{p) = 1 

and ?2('^) ~ '2('^ + tt) — 2 sin | cos 1^, with (p being an 
angle between y and m. We then have: 



m • y 



= -d-TTsin -. 



Q r-2TT l-hiv) 

-d^y = — 2 sin — / d(p 

2 Jq Jo 



dy 

cos ip 

(B6) 



For the integral I3 we have: 



1 



1 



- 2s ■ y 
m • y 



2 sin I Jo cos (p Jo 

^ I o 9 \ . •) 9 
= -2tt cos-' - + - sm - 



dy 



y + cos — — 2y cos - sm ip 



(B7) 
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where we used the relation hi^) = hi—'f)- 

Using integrals h dBS]), 12 die]), and h (EZl in 
Eq. (|Bip . we obtain equation ([51)1 : 



37r + 27rsin 



61 /4 



2 V3 



■ — In I tan ■ 



Turning to the variable z = — sin^(6'/2) the integral I[ 
is reduced to 

^ Jo J- sin2(e/2) ^ 

with 



Appendix C: Calculation of the interaction function 



f: 



(2) 



" (2) 

The interaction function F2 is the second variational 
derivative of the many-body contribution to the intcrac- 

~ ^2^ 

tion energy, E2 (I83p . with respect to the momentum 
distribution. It reads: 



(k,k') 



2r 



where 
/(=2/ 

J|ki 



d^fci |k-ki| ■ |k^-ki| 

k2+k'2-k2-kl 

d^fci |k-k'| • |ki-k' 



k,|<fe, k2+k?-k'2-k2 



(^k-fJc'-ki-ka I 
<5k 



2 0k-*i-k'-k2 



(CI) 

(C2) 
(C3) 



The integration area for /( is shown in Fig.[5J where the 
distance between the points Oi and P is RoiP = k/2fci^, 
RpOa = k'/2fci?, RoiTV = ki/2fci?, and Row = x/2. 
We thus have Rnp = (k — ki)/2A;F and Rnp' — (k' — 
ki)/2fci;-. In the region of integration we should have 
\Ro,n\ = ki/2kF < 1/2. This leads to 



IRwpI ■ IRwp' I ,2 



2tt rhiv) 

d'^'n = — / dip I xdx 

"'0 



sin2(6l/2)]2 - 4a;2 sin2(6l/2) cos^ ip 



\e/2) 



(C4) 



where ip> is the angle between x and m (see Fig. [2]), and 
the quantity ^3(1^9) obeys the equation 



^3 ( v) ~ 2 cos - sin (p ■ I3 {(p) + cos^ - 



i?, = + 4z sin^ - sin^ ip + 4 sin^ - sin^ tp. 




FIG. 2. (color online). Left: The integration area for I[ (in 
blue): RoiP ~ k/2fcF, RPO2 ~ k'/2fcF, RoiJV = ki/2fcF, 
and if is the angle between the vectors Rop and Ron, which 
is the same as the angle between m and x. Right: The inte- 
gration area for I2 (in red): RqiP = k/2fcF, R02P = k'/2fcF, 
Ro2iv = ki/2fcF, a is the angle between Rpm and Rpjv, and 
<j!> is the angle between Rpjv and Roo2- 



It is easy to see that: 



rll{v)-sin^{e/2) ^ 
Ir — I dz 

-sin2(e/2) ^ 



: |V^- y/aln (^2a + bz + 2VaR 
^ In (2VR + 2z -f I 



i^(v)-sin^(9/2) 
-sin2(6l/2) 



i^(v)-sin2(e/2) 
sin2(e/2) 



— Irt — Iri, 



where a — 4 sin^(0/2) sin^ b = 4sin2(6'/2) sin^ and 
the symbol P stands for the principal value of the in- 
tegral. The quantities Ir^ and Iri denote the values of 
the integral at the upper and lower bounds, respectively 
(in the last line we have to take the principal value of 
the integral and, hence, if the upper bound of the inte- 
gral is positive we have to replace the lower bound with 
sin^(6'/2)). Then 1,.^ and /,.^ are given by: 
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Ir^ — 2| siniy9| • l{(p) — 2 sin^ — | siniy9| 



In ( 8 sin^ - sin^ ipj + In f sin^ - + cos - sin ip ■ l{ip) + l{ip) 



+ 2 sin^ -| siniy9| • In |2 cos - sirn^ • l{(p) \ + 2 sin'^ - sin^ (p 



Iri = sin^ 2 ~ ^ ^"^^^ 2 ' 



+ 2 sin — sin (/s • In ( 4 sin — sin ip 



ln4 + In ( I siniysj • l{Lp) + cos — siniy9 • l{Lp) + sin'^ — sin^ ip 



In I 4 sin"' - J + In (sin^ ip + | sin(^|) 



.n(7.. . I ■ 7 ^ 

2 sm - 1 sm (/?! • In I sm - 



r 



The integral /( can be expressed as: 



1 

2 Jo 



[Ir-l - Irl]dip, 



l^{ip + tt) = sin^ |, h{ip) — h{(p + tt) = 2 cos | simp, and 



h{^p) + + tt) = 2y/cos2 I sin^ + sin' 
obtain: 



2 fl 



Wc then 



and for performing the calculations we notice that /3((^) • 



= - sin^ - 7rln2 + 7r/2 - tt In sin - + 41n | cos -| - 41n(l + sin-) +g{6) 



2ti 



4 arcsin(sin |) 



TT — 2 arcsin(sin — ) + | sin ( 



(C6) 



with 



— I 2 sin^ ip In I sin f + ^ sin^ - + cos^ - sin^ ip j dip. 



(C7) 



The integration area for 1 2 is shown in Fig. [21 and we 
get: 



4 



|m| • IRpATi 2 



m • y 



cos 6 



(C8) 



where we denote Rpjv — P, and (j) — a — 9/2 is the angle 
between the vectors m and p, with a being the angle 
between the vectors RpM and Rpat (see Fig. We 



then have: 



TV /'Sin a 



pdp 



!'. = -% da , 

2 ' /o cos(a-6'/2) 

2^*, l + sin(6l/2) ^ . 6'' 

cos^ - In + 2 sm - 

2 1 - sin(6'/2 2 



(C9) 



Using I[ den and ^ (IC9| in Eq. jCT]) wc get equation 

^ ('2) 

5]) for the interaction function F2 (0). 
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